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The formulae of particle current as well as spin- and angular momentum currents 
are studied for most spin-orbit coupling (SOC) systems. It is shown that the con- 
ventional expression of currents in some literatures are not complete for some SOC 
systems. The particle current in Dresselhaus system must have extra terms in ad- 
ditional to the conventional one, but no extra term for Rashba, Luttinger model. 
Further more, we also prove that the extra terms of total angular momentum appear 
in Rashba current in addition to conventional one. 

I. INTRODUCTION 

In quantum mechanics calculating current is crucial to applications. A widely accepted 
approach is using the correspondence regulation from classical to quantum mechanics. For 
example, the formula of charge current density is introduced as 
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t) =ei?e{^p„^J,(r,t)^^„(r,t)}, 



where e is an electron (or hole) charge, 1^ the velocity operator defined by ^ = j^pr^^H], 
H is the Hamiltonian, \l/„(^,t) is the wave function of an admissible state of the 
Hamiltonian system and p„ is the probability of \I'„(^,t) appeared in the quantum mixed 
state. For example, if the system is in equilibrium, the distribution of the probability is 
Pn = e~^^", where /3 = l/(fcT), En the energy of state \l/„(T^,t). Summation n is taken 
to all admissible states of Hamiltonian. Further, the definition of charge current density is 
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extended to define some other currents like spin current density 

n 

and orbital angular momentum current density 




where P = x "p*)^ is orbital angular momentum along z direction. And symmetrization 
in the order of ^ and P (or P) has been apphed. The above conventional formula have 
been widely accepted as natural definition of currents in literatures, especially in the ones 
of semiconductor systems including spin-orbit-coupling (SOC). 

In this paper, we will prove that the conventional formula of particle current are not 
generally correct. They can not satisfy the continuity equation for some systems. Based 
on the continuity equation, [9] we start our calculation from dwA^J^ ■,t)/dt for observable 
quantity A, as will be illustrated in Sec.2, and utilize the continuity equation to obtain 
the definition of ^-current j a- We can find in the case hke 2d Rashba, 3d Luttinger and 
the spin independent Hamiltonian H = 2m* ~^ ^(^)) where y(T^) can be an arbitrary 
position dependent potential, the conventional particle current formula are still corrected. 
However, we prove that the particle current will have a non-trivial additional term except 
conventional one in some semiconductor, such as the system with Dresselhaus SOC and other 
model Hamiltonians including the terms p*^" of momentum operator with power n > 2. We 
will derive the expressions of the extra term for these specific systems. 

This paper will focus on some systems with spin-orbit-coupling such as 2d Rashba[l], 
3d Luttinger [2] as well as Dresselhaus systems [3] that attracted great interest since the 
experimental demonstration [4-6] of spin hall effect in some semiconductor materials. Our 
study in this letter is based on non-relativistic quantum mechanics. 



3 



II. PARTICLE CURRENT FOR SOME COMPLEX SYSTEMS 

In general, we study an observable quantity A. Prom quantum mechanics, A corresponds 
to a Hermitian operator A, then 

(A) = Tr{lp} = rr^{l|*„)p„(*J} 

n 

n n 

where we have used (^|y4|^') = A{^')5{1^ — ~r') which is diagonal in position representa- 
tion for the most of cases in quantum mechanics, A(l^') is the operator A in position repre- 
sentation and simply expressed by A. p{= n) Pni"^ n\) is the density matrix of a mixed 

state, Pn (g [0, 1]) the probability of appeared in the mixed state and = 1. \^ri) 

n 

is assumed here to satisfy the time dependent Schrodinger equation ihd{\'^ n)) / dt — 
Since {A) is real, so the above equation is equivalent to 

{A) = ReJ2pn dV{<lfn\T'){T'\A\<lfn) = / WA{T',t)dV (1) 

and 

WAiy.t) = ^PnRe{'^n\'r){'r\A\^n) 

n 

= Y,PnRe{^i{V,t)A^n{y,t)], 
n 

n 

Where WAiy it) is defined as the density of quantity (^4) at time t and position 7^. When 
[1,^] = 0, 

n 

= ^7^PnRe{'^n\AH-EA\^n)=^ 
n 
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A is a conserved quantity, therefore it yields a continuity equation: 

dwA{1^,t) -> 
dt =-^-^^- 

Particle number is one of conserved quantity, here it is relevant to choose A — NI, I is an 
identity operator. In general, one define a conserved current related to mechanical quantity 
A from above continuity equation. Now the key point is to find out the calculation formula 
of current j a ^or each specific Hamiltonian. When A is relevant to the total particle number 
N,A — NI, I is identity operator, it yields a definition of particle current j . Total number 
of charges eN is also a conserved quantity, the relevant operator A — eNI and conserved 
charge current j e — ^ j ■ In general, A could be others, so one can define other currents. 
In this note, we only consider mixed state with time independent probability p„. A is also 
studied within time independent. 

For the particle with spin s = 1/2, the state |^') can be described by 

l^i) 





tpii r ,t) 

where ipliy ,t) is complex conjugation of ipi{'r^,t) {i = 1,2). The Schrodinger (or Pauli ) 
equation is 

where Hamiltonian is a 2 x 2 matrix. Since now, for simplifying notation, we will use ^ 
to denote '^{'r',t), \E'^ to \E''I"(T*, t). The time evolution of quantity w^^\'r',t) is 

In equation (2), the operator A has been represented in position space. Meanwhile, it must 
bear in mind that and are the time and position dependent wave function. When 
A — I2 (a2x2 identity matrix), 

n 
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nCj^jt) is the particle density, so the following continuity equation defines the particle 
current density j . 

It describes the conservation of particle number. For a time independent Hamiltonian, 
the wave function = exp(-i£;„i/fi)*„(0), -H^*n(0) = -E„*„(0). Therefore we have 

dn{r, t)/dt = 0. It reduces to V • ^(y, t) = 0. 

For many Hamiltonians without spin-orbit-coupling, H = Hq + ¥{1^), Hq = /2m* — 
— (/i^/2m*) V^. In terms of the equality Hq = — |V ■ ["r*, ifg] = ' Ij^iH]^ equation (1) 

can be easily calculated as 

dt \v^;i / "in ° " 

= -V • i?e ^o]^n} = -V • i?e ^]^n} 



(3) 



In equation (3), we have used the real property of potential ^(^) that results 
+ ^l^-rj:y{~r)'i/n — 0. Thus we obtain the standard continuity equation 
and the formula of particle current 



dn{l^ , t) 
dt 



(4) 



Equation (4) is called as " conventional formula" . The formula is general for any potential 
Vi^r) in Hamiltonian H if the potential is position dependent only. However, if potential 
V contains momentum operator like the case appeared in some effective Hamiltonians of 
semiconductors, it will be different and the above current formula may need to be confirmed 
or to be modified to including some extra terms. In recent literatures, ones usually extend 
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it to define other kinds of conventional current like spin current, A — S 



2 ^ ' 



(5) 



, but not conserved, where symmetric presentation is applied for two non-commutating 
operators ^ and S . 

Now we discuss the particle current in some complex systems with SOC in which the 
momentum operator appears in "potential" of Hamiltonian. First, we study the current 
formula for the Dresselhaus Hamiltonian that contains triple power of momentum operator. 
The Hamiltonian is 



-f^o + He 
1 



H 
Ho 

Hn = V\P. {pI - pI) a. + Py {pI - pI) + p. (pI - pi) a,], 



2m* 



p^ + V{r 



where {(Tj : i = x,y,z} are Pauli matrices. The question is whether the particle current can 
still be calculated by conventional formula (4). We rewrite the Hamiltonian as 



H^Ho + Hd = -- V • [r. Ho] + ViV) + Hd, 



It yields 



dn{^ , t) 
di 



Re^pn 

n 

Re^Pn 

n 

+Re ^ pn 



ih 



Similar to the deduction of equation (3) for the terms having Ho in right side of above 
equation, we directly have 



dn{/r , t) 
dt 



= -V-ReY,Pn 



+Re pn 



1 



^li-HD^nj + ij^HD^n] 



t 



(6) 
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If the second term in right side of above equation equals to — W- Re Pn "^ni j^ir, -f^oj^n ) 
the final formula of particle current is just the conventional one. But, in deed, they are 
different. Appendix (A) gives the detail proof to show that an extra term must exist in 
addition to the conventional current formula. The continuity equation should be 

dn{~r , t) 



y ' \ J conv + J extra): 



dt 



j extra ^ ^ Pn j extra i"^) 

n 

G^tra)x = \h'v{{^l-Vl)[^i{a^^n)] 

-2V.V,[*t (<7,*J] + 2V.V,[*t (<7,*„)])} 
G^tra)y = l^Mi'^l-'^lmU^y'^n)] 

-2V.V.[*t C^x^n)] + 2V.V,[*], {cTy^n)])} 

where j conv is the conventional formula of particle current density. In appendix A, we also 
prove that the term j extra is real and position dependent if the position dependent potential 
VCr^) is included in Hamiltonian. Therefore, V ■ j extra 7^ 0, so the nontrivial extra term of 
particle current do appear in Dresselhaus system. 

We also checked the formulae of particle current for 2d Rashba 
(h ^ Ho + aZ ■ {'p X ~a)^ and 3d Luttinger [h = HqU + K'P ' ^)^) model Hamilto- 
nians, where 5" is the 3/2 angular momentum operator and /4 is a 4 x 4 identity matrix. 
It is found that no extra term of particle current appears for these two models. The main 
source to induce an extra term of particle current for Dresselhaus Hamiltonian seems from 
the triple power of momentum operator in its Hamiltonian. Indeed, in general, if the 
Hamiltonian includes momentum operator with order higher than power 2, the extra term 
of particle current will appear. For Rashba and Luttinger models, the highest power of 
momentum operator is 2, but Dresselhaus is 3. So this is the reason why the extra term 
of particle current appears in Dresselhaus, but not in Rashba and Luttinger systems. In 
section 5, we will show that from an extended Noether's theorem. 



8 



For showing such argument more clearly, we study a simple model Hamiltonian without 
SOC: 

H = Ho + Hi + V{r), 

Hi = Pp", (8) 

where /9 is a coupling constant. In section 4, we will show that from the extended Noether's 
thoerem 

dn{l^ ,t) ,— > — > , 

77, ^ ' \ 3 conv ~l~ 3 extra ) i 

at 

y extra = 2ph' ^ p„/m{ V (^J^V^^n) } (9) 

n 

For nonhomogeneous system, \E']jV^\l/n is complex including both real and imaginary parts in 
general, j conv and j extra in equation (9) are position dependent, V • j conv and V • j extra 
can be nonzero in general case. Therefore, they do not satisfy the continuity equation 
separately. But the combination of them gives the following continuity equation for any 
stead state (^^^ = 0). 

V • = V • ^conv + V • y extra = (10) 

In this example, we clearly show that the extra term of particle current comes from the part 
of Hamiltonian with momentum (n = 4) . 

III. SPIN- AND TOTAL ANGULAR MOMENTUM CURRENT 

When the operator A in expression (1) is spin or the total angular momentum (TAM), 
— P + S^, equation (2) defines the time evolution of spin-density or the TAM density 
that could relate to the spin-current density or total angular momentum current density 
(TAMCD) both often appeared in recent literatures [7-10]. For the Rashba Hamiltonian, 
[S^, Hr]- ^ 0, the spin is not a conserved quantity. Ref.[8,9] have shown that 

dt dt^ 

= -V • 7p + Tir^ (11) 



9 



j s — ^ Pn 3 is so called "spin current" widely in recent literatures, 

n 

lu. = 5^Pn7L"^ = Re I J^Pn^'I ((«(? X ^) X ^)^r. 
n \ n 

which is called as spin "torque". The conventional spin current density j s just partly 
contribute to the time evolution of spin density 10^^(1^ ,t){— ^ pnw'^^^(/r ^1)). Even 

dw^^Cr^ ,t)/dt = for a steady state, the conventional spin current density j s is not 
conserved ( V • j conv 7^ ) and satisfies the equation V • j conv = j u- 

Now we study the current of TAM = + in Rashba system, where is the z 
component of orbital angular momentum. The case is different from particle current. In 
case of particle current, the number of particles is conserved, but relative operator of particle 
number is momentum operator independent, A — NI. Here is also a conserved quantity 
since [J^,if]_ = for Rashba Hamiltonian, but it contains the momentum operator. So 
we would check whether its related current has extra term even for Rashba Hamiltonian 
without higher than 2 power of momentum operator. After long careful calculations (see 
appendix C), we obtain 

dtY,Pnwf^{r,t)=dtY,PnRe[^ir^n]^-'^-lL, 



where 



J L 



1 



1 

n 

(^v, - yV,) V^n + {{xVy - yV,) V^n^^n} (12) 
The conventional definition for current density in literatures is defined by 

K L J J 

After taking the real part, it becomes: 

7conv = lReY,pn{'^i ^[r,H]r + r^[r,H] 
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Now we use the following relation 

= ih[{py/m* - aa^Ye ^ - {px/m* + aay)~ey] 

then change the form of j conv into 

n 

+]^Re^pn [iMiiipy/m* - aax)~ex - {Px/m* + aay)'ey]^n] 

n 

= i?e5]p„{^l^[^,^]^^^n} + ^i?e5^p4^^*l[^^.-^^.]^n} (13) 

n n 

In the last step of above equation, we have used that ^iJa^^n are real since {a^, ay, a^} are 
Hermitian and trace has been taken in spin space ( but be care for that no trace taken in 
position space). In appendix C, we obtain 

7 extra ^l^L- 7 conv ^-^Yl ^"^'^(^n*") ^ ^ (1^) 

n 

Thus, we conclude that conventional formula of TAMCD, j conv, is not complete. It must 
have an extra term even for Rashba system. Therefore, we suggest that one should be careful 
in using the conventional current formula and do the deduction of formula carefully for every 
new type Hamiltonian, particularly for some Hamiltonians with SOC in semiconductors. In 
this note, the TAM, (Jf), is not a conserved quantity for Drcsselhaus Hamiltonian since 
HDresseihaus] - 7^ 0. So, wc did not attempt to check the existence of its extra part. 

IV. NOETHER'S THEOREM AND EXTRA TERM OF CURRENT 

Based on time dependent Schrodinger equations and the particle number conservation 
(see section 2), we have derived new expressions of conserved particle current for some 
Hamiltonians, which include term that the order of momentum operators in it is higher 
than 2. In parallel to that, in fact, expressions of particle current can also be derived by 
Noether's theorem in field theory from U (1) gauge invariance of Lagrangian. However, 
the usual Noether's theorem only consider systems with term /2m in the momentum 
dependent part of Hamiltonian. So the Lagrangian of the system is only a functional of field 
functions and their first order derivatives, 0(x^), (f)\x^), d^(f){xij) and 5^0''' (x^), for a complex 
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scalar field, where x^j, — {t, x,y,z) , fj, — 0, 1, 2, 3; {do, di, 82, 83) — {df, d^, dy, d^). Simply, we 
denote this Lagrange as £(0, 0^, 9^0, df^cf)^). Let us briefly recall the usual Noether theorem. 
An action is defined by 

S ^ I dx'^£{(f), (f)\ 9^0, d^(f)^), dx'^ = JJ^ dx^,. 



(15) 



Where two variational functions 50 and 50^ are independent. According to the least action 
principle, 5S — 0, which gives two conjugate dynamical equations of and 0^ One is from 
the term of 50 and another from 50^ For a complex scalar field, the Lagrangian is usually 
given by 

£(0,0t,a^0,a^0t) = 0tzao0-2^ (^,0)^^,0) - 0tv({a:^})0, 

i=l,2,3 

where h—1. From equation (15), we have 

dC ( ^ dC 

- o„. 



2m 

It is Schrodinger equation: 



a0t \^''d{d^cj>^)^ 

ido<j> + ^d^<j>-V{{x,})<f> = 0. 



^|0 = H<l>,^ 

H^-^^Vi{x,}l 



Equation (15) also yields an equation: 



E d,F,^Q. (16) 

M=0,l,2,3 



which defines a AD current: 
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Input usual Lagrangian L into equation (17), we have 



Fi = -2^P,0)^(5(/) + (a,(/))5(/)t],i = 1,2,3. 

By t/(l) gauge transformation 

e'"0 (1 + iQ;)0 = + iQ;0 ^ 
= iacj), 6(j)'^ — —ia(f)\ 

and due to U{1) gauge invariance of the Hamiltonian, we obtain 

Fi = ~^{{di<t>y<P-<P^di<f>}^-Re<f>hdi<f> 
2m ^ ' m 



—aRe ^4>^-[ri, H]_4>^ = —aRe {^^Vj^j , i = 1, 2, 3. 



We obtain the continuity equation of particle 



dt 

This is just the conventional expression of particle current. 



+ V-i?e|0^^[r',ii']_0| =0. (18) 



Re{c^^-[1^,H]_<P\ (19) 



For Rashba and Luttinger systems, the usual Noether's theorem can be applied since the 
additional spin-orbit-coupling term in these two Hamiltonians are momentum operator p 
or dependent. Therefore, it can be unstood that why their expressions of the particle 
currents are also conventional. It is consistent with our conclusion. However, when we study 
the system like H = p^ /2m-\-j3p^ or Horesseihaus, the terms of high order derivatives, inherited 
from Hamiltonian, should also be included in the Lagrangian: £(</>, 9^0, d^d^cp, (0 0^))) 
thus the usual Noether's theorem should be extended. According to variational method, 

SS ^ [ dx"5£([0,a^0, 9^9,0, (0^0^)]) 
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where dx^ — Y[ d^n- Now the 4D current is defined by 



+ --- + {4>^4>') (20) 



and 



9£ ^ f dC \ ^ ^ ( dC 



A£ = — - 5„ ( ,^7.^ 1 + d,d^ 



d(t> '\d{d^<l>)J ^\d{d^dA) 



Thus 



5S^ I dx^AC5(t)+ I dx'^d^F^^O. 

JVi Jv4 

According to the least action principle, a suitable definition of Lagrangian is necessary for 
obtaining the correct dynamical equation of from AC = and conservation law from 
dtj,Fn — 0. Therefore, the expression of particle current j can be obtained from the 
invariance of U{1) gauge transformation for such extended system: 

+ ... + (4,^4,1) (21) 
Comparing equation (21) with (15), there are some extra terms F^^*''" : 

+ . . . + (0 ^ 0t) 

This is the reason why extra terms do exist in particle current for the model and p^- 
Dresselhaus Hamiltonian. In this paper, as an example, we use the above extended Noether's 
theorem to study the p'^ model system in detail only. Further study on Noether theorem for 
other systems must be interesting and it is ongoing. 

The central point of Noether's theorem is to find out a suitable functional expression of 
Lagrangian, then one can derive the time dependent Schrodinger equation and conservation 
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law. For mode, we choose the following symmetric Lagrangian 

£[0 (x) , d.cP {x) , dl4> {x) , {x) , dldlcP {x) ; (0 ^ 0t)] 

-(a,0)^(a,aj0) + 0t(a2aj0)], (22) 

After long algebra, but no difficulty, we obtain Schrodinger equation 



dt \2m 

and the following equation for AD conserved current 



dC dC dC 

= ^am'"^ * am' " '""am*'* 

-5{d.d,<t>)-{d,-j^^)8{dA) 



d{d,dl<f) ^ ^ ^ 'd{d,dlcf) 

dC dC 

and conservation law d^F^ = 0. It is an extension of the conventional expression = 
-g^^S(j)+h.c. for usual Hamiltonian. Now applying U (1) gauge symmetry of the Lagrangian 
to eq.(22), performing transformation: — > (j) + iacf), (j)^ — > 0^ — ia(f)\ where a is a small 
constant, we gain the expression of conserved particle current 

Icons = (V0) - (V0)^ 0] + z2/?[0t - (VV^^t) 0] 

-i [-5 (VVV^) + (VV^ 

+3 (V0t) (VV) + 2 {V(t>f ■ (VV0) 

-3 (V V) ^ ( V0) - 2 (VV(/.^) • ( V0)] (24) 
which is essentially the same as eq.(9) in Section 2. We have proved that V • 

^ 3 cons ( J conv ~l~ j extra)^ 0- 
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We also have proved that a suitable symmetric Lagrangian for Dresselhaus model 
can be found and correct Schrodinger equation and particle current are gained. The ex- 
pression of particle current is essentially the same as the one in section 2 with a difference 

A i Dresselhaus = V X A , 8. TOt VCCtor, Only, SO V • (^A j Dresselhaus^ = 0. Fuithci StudicS On 

the extension of Noether's theorem about the construction of Lagrangian for such unusual 
system, the multiple correspondence from Hamiltonian to Lagrangian and finding general 
approach is our next studies. 

In summary: we have studied the conventional formula of particle density for some Hamil- 
tonians such as Rashba, Luttinger, Dresselhaus Hamiltonian as well as a model Hamil- 
tonian. It is shown that the conventional formula of particle current density is correct for 
Rashba and Luttinger Hamiltonians. But it lacks an additional term of current density 
j extra foi p'^ Drcssclhaus Hamiltonian and p'^ one without SOC. The spin and TAM cur- 
rents arc also addressed for Rashba Hamiltonian. Some detailed proofs arc presented in 
appendices to show the existence of extra terms of currents. The results point out that the 
conventional expression of TAMCD is not complete even for Rashba term. It against the 
conservation law of the total angular momentum if the extra term is neglected. In deed, 
the extra term saves the conservation of total angular momentum for Rashba system. In 
the final, a brief explanation of the extension of Noether's theorem is presented. It just 
shows that the deviation of the expression of particle current from conventional formula 
derived by Noether's theorem for some system with SOC is reasonable and no doubt. The 
general extension of Noether's theorem is still not quite clear and in studying. However, 
our approach to derive the conserved currents is based on Hamiltonian and time dependent 
Schrodinger equation, our results are rigorous and unique up to a rot vector (divergence =0) 
without any doubt. A correct formula of particle current is important in studying transport 
problems in recent SOC systems. 

Acknowledgement: The work is supported by National Natural Science Foundation of 
China (No. 10674027) and 973 Project of China (No. 2002CB613504). 
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V. APPENDIX 

A. Extra term of particle current for Dresselhaus system 

In this part, we would show some details leading to the extra term in the current density 
for Dresselhaus Hamiltonian H — Hq + Hd- Let's begin from the equation of (6) in the 
body of the text 

r9T7 1 1 ^ 11 ^ 

at = -<2^^ ■ - . [T>, /f„l)*„) 

+ (^fiD*«(7', *)) t„(T,t) + (^^SBt„(7', t) j . (A.l) 

Now we focus on the deduction of the last two terms of Eq. (A.l). Considering SOC part 
of Hamiltonian 

Hd = V [Px {pI - pI) cTx + Py {pI - pI) CTy + Pz {pI - pI) (Tz\ , 



we obtain 



1 

ih 



x,Hd 



+ 2r; {pyPx(Jy - PzPx<^z) ) =v{pI- pI) (^x, 



(A.2) 



1 

ih 



y,H 



D 



+ 2r/ {PzPyCTz - PxPy(^x) = V {pI - pI) <^y, 



(A.3) 



1 

ih 



D 



+ 277 {PxPzCTx - PyPzCTy) ) ^ V {pI - P^y) 



(A.4) 



H 



D 



ih 



+ '^V{Py(^y-Pz(yz)pl 



+277 iPz(7z - Px(7x) pI + 277 (PxO-x - Py(7y) P, 



P 



ih 



H 



D 



-I- 



r,Hn 
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Then we have 



1 

3 



4v 



+- 



ih 



V- 

1 

ih 



t 



1 



t 



1 



1 

ih 



r,HD 
1 

ih 



D 



Inserting above result into the equation (A.l), wc obtain 



dn 
'dt 



1 



Hh 



1 



A 7 



+ -Re \ (V* 



1 

ih 



Now we try to express |i?e 



by V • jarfd; then the analytic 



formula of extra contribution to particle current can be obtained. By using (A. 2), (A. 2) and 
(A. 4), j^I't'^Hd] can be substituted as following 



1 

3 



ih 



ih 



-l^^Miii^l - V') (Jy - 2V,Vya, + 2V,V^(7,] Vy^n} 

-h''7]{Vy^i [{Vl - Vl) ay - 2V,V,(7, + 2V,Vya,] ^^j 
-h''7]{V,^i [{Vl - Vl) a, - 2V,V,c7, + 2VyV,ay] 



It can be changed to 



1 
3 



+ (2V^V2,(7^*n)^ Vj,*„ + {Vy-^i)2V^Vya^'^n 

- (2V,V,(7,*„)^ V,*„ - (V,*t )2V,V,(7,* J + CP. 
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Where CP. is cyclic permutation of index {x, y, z}. Above equation can be changed to 

-V,[(V,(7,*„)W,*„]} + aP. 



-Vx[{Vx(yy^n)^Vy^n + (Vj,*],) (V^(7j,*„)] 
+V,[(V,a,*„)^ V,*„ + (V,*y (V,(7,*„)]} + CP. 

where 
So 

7^(n) _ ~^{n) , ^/(n) , "^(n) _ ^(n) , ^(n) /a 
J ~ J coni; ^ J ^ —J conv ^ J extra l^-^J 

where 

J extra J ^ 
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And 



= -^{^n(^ {{pI - pI) - '^PyPx(Ty + '^PzPx(^z) 

+ (^ {{pI-pV) - '^PyPx<7y + '^PzPx<7z) *n)'^*n} 

So the X component of the extra particle current density is 

I J extrajx — -^^x ^ J x 

- 2V,V,[*t (^^^^)] + 2V,V,[*t (^^^^)]| 
1 



^'^{(V^ - V^)[*t (<7.*„)] - 2V.V,[*t (<7,*„)] 



6 

+ 2V,V,[^t (^.^n)])} (A.6) 
And we can obtain other components of j ^extra by cychc permutation of indices. 

- 2VyV,[^Ua,^n)] + 2VyVM i^x^n)])} (A.7) 



- 2V.V.[*t (<7,*„)] + 2V.V,[*t (<7.*n)])} (A.8) 

which is the extra term of the current density j extra in Dresselhaus Hamihonian as stated in 
the body of the text. It is not zero in general since is real and position dependent 

for general Hamiltonian. Therefore the terms in above equations are real and not zero after 
the derivation of position in general. 
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B. Extra term of particle current for Hj = (3p^ system 
H = Ho + Hi + V{r), 

And since 

= /5 {(v^t) . ^3^^ ^ (i?^„(r))t . (v^„)} 

= {(V*t) • (V3*„) - (V3*„)^ . (V*„)} 

= V- {(v*t) (v^*„) - (v^*„)^(v*„)}, 

we have 

= -V-it:e{*t^[r,^o]*n} 

+v ■p{-ihf[{v^i) (V'^^n) - (V'*„)^ (V^.)} 
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So we get the conserved current 



t = i?e{*tl[r,#o]*n} 
in 



Compare with the conventional definition, 

~Jconv = Re{'ifi^[y,Ho]^n} 

we find the extra term of particle current appears. 

j extra j j conv r>"^^{^n"^[ ^ ' ^l\^n} 

2 in 

= -/3(-^n)^{vi/|^vv2M/„ + (vvi/y (v^^n) 

= 2f3K'lm{W {^^V^^n)} 
And in general case, the state can be expressed as 

7 extra = -/3in^{*J,VV2*„ + (V*],) (V^*^) - (V'*„)^ (V*„) - (VV'*„)^*„} 

= -/^^^' E E ^*(^- *)^(^- 0[^^n^n - ^^^^ - iK'~k^ + i^;2^]e-(^"-^)- 
^ 

We reached our conclusion that, in general case, includes many components 

{A{k n,t)e'^'^'^''^ }, therefore the extra charge current density j extra is position dependent 
and not zero. 
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C. Extra term of total angular momentum current for Rashba system 

First, we have Rashba Hamiltonian: H — Hq + Hr, and it is easy to check following 
relations: 



ih 



-H, 







ih ^ 



1 ^ 
— p 

m* 



CXCTy 6 X ^ y 



■V.|l[r,5„: 



The time evolution of the density of TAM satisfies following equation 



where, in the second step, we used the fact that commutes with the Hamiltonian H. 



-V 



- V- 



ih 



V J^^, 



Considering 



h 



1 / h 



2 V irri^ 
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t 



= V • (^(aay-^, - aa^y) >*n) } - V • |^_L[7>, | 
Where we have used the following relation 

Therefore, we obtain 



a(*t (-?,*) J'*„)/a( = -iv ■ ('*' 111*, ffoi - iv 



-V- 



t 



Since ^[7", H'r] = aay~e^ - aa^~ey, it yields 



Thus, we have 



-V- 



Meanwhile 



— p + auy e X - ctcTx e y H'r 
m* ' 



h 



im 

h 



V + auy e X - aa^ e H/„ 



im* 
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It results 



8(*iJ'*„)/Si = 



n V n J 



And use the fact that ^l^ai^'l is real, we have 

r ^ 



-Re 



2im 



1 h'' 



1 /i^ 

y ^y^ X 



Therefore, 



1 /i- 



+ 



2 2m 
2 2m' 



+ {xVy - yVx) V*„ + {{xVy - yVx) V*„)^*„} 



The conventional current of TAM J£: 



(C.2) 



(C.3) 



ii?e5]p„*t |^[- ^]> + j._L[7^,^]U^ (C.4) 
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Due to the following relation for Rashba Hamiltonian: 



ih[{py/m* - aa:c)~e^ - {px/m* + aay)~ey] 



(C.5) 



we have 



+ ^^e^pn {im\\{pylm*)'e ^ - (pa;/m*)'?y]*„} 

n 

1 

+ ^— {*UVy"?x - V^"?2/)*n + ((V^,"?^ - V^"?j,)*„)^ (C.6) 

where, in the second step, the property that '^\(Ji^\ is real is again used. Compare (C.6) 
with (C.3), we have 

J extra J L J conv 

1 

So the components of j ''^]tra 

1 



■ {n)x 
J extra 



2 2m* 



1 



4m^ 



Similarly, we have 



(C.7) 
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So finally we arrive at 

~t extra = ^ Pnl'^ {^1^ r^) (C.8) 

n 

Where ^l^^n is the density of particles and real. It must be position dependent in general, 

thus 7extra = £*[xV y " yV,.]V(*t*n) ^ 0. 
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